hep-th / yymmnn 
March 2008 



One-loop corrections to the instanton 
transition in the Abelian Higgs model: 
Gel'fand-Yaglom and Green's function 

methods 



Jiirgen Baackc 



Fachbereich Physik, Technische Universitat Dortmund 
D - 44221 Dortmund, Germany 



1 e-mail : j uergen . baacke@ tu-dort mund . dc 



Abstract 



The fluctuation determinant, the preexponential factor for the instanton 
transition, has been computed several years ago in the Abelian Higgs model, 
using a method based on integrating the Euclidean Green' function. A more 
elegant method for computing functional determinants, using the Gel'fand- 
Yaglom theorem, has been applied recently to a variety of systems. This 
method runs into difficulties if the background field has nontrivial topology, 
as is the case for the instanton in the Abelian Higgs model. A shift in thre 
effective centrifugal barriers makes the s-wave contribution infinite, an infin- 
ity that is compensated by the summation over the other partial waves. This 
requires some modifications of the Gel'fand-Yaglom method which are the 
main subject of this work. We present here both, the Green' s function and 
the Gel'fand-Yaglom method and compare the numerical results in detail. 
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1 Introduction 



The computation of functional determinants for various background field 
configurations has recently found renewed interest. In particular, the ele- 
gant Gel'fand-Yaglom [H El EJ H] approach, denoted sometimes as "the Cole- 
man method" (as presented in [5]), has been considered by various authors 
[SI El El [TOl [HI [12], following earlier work on vacuum decay and bub- 
ble nucleation [131 O EES! EES1 EH EE]- Another method for computing the 
functional determinant, based on the integration of the Euclidean Green's 
function, has been used in Refs. [19] and [201 EI] for computing the fluc- 
tuation determinants for the instanton in the Abelian Higgs model in 1 + 1 
dimensions and for the sphaleron in the SU(2) Higgs model in 3 dimensions, 
respectively. When one tries to naively apply the Gel'fand-Yaglom method 
to these cases of topological solutions in gauge theories one encounters a spe- 
cific difficulty: the topological soliton modifies the centrifugal barriers and as 
a consequence the contribution of the s-wave sector is infinite. The problem 
is avoided when using the Green's function method if the summation over 
partial waves is carried out before the integration of the Green's function. 
For the Gel'fand-Yaglom method a modified approach is required; it is the 
aim of this work to elucidate the problem and its solution. At the same time 
we consider the relation between the methods, analytically and numerically. 
We use the Abelian Higgs model here mainly as a typical example of a model 
with a topological soliton, the problem is present for the sphaleron transition 
in the electroweak SU(2) Higgs model as well. 

The Abelian Higgs model in (1 + 1) has found considerable attention since 
on the one hand it shares certain features with the electroweak theory, and 
on the other hand it is simple enough to serve as a theoretical and numerical 
laboratory. In the context of the baryon number violation the high temper- 
ature sphaleron transition in this model has been studied [22"l [23] for which 
exact classical solutions and an exact expression of the sphaleron determinant 
[2"H [25] are known, thus providing a complete one-loop semiclassical transi- 
tion rate which can be studied numerically on the lattice, e.g. by measuring 
the fluctuations of the Chern-Simons number. 

Another prominent feature of the model is the existence of instanton so- 
lutions [261 EI] which give rise again to fluctuations in the topological charge 
of the vacuum and thereby to baryon number violation. In the dilute gas 
approximation for the instantons transition rate, or equivalently the density 
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of instantons in the Euclidean plane, is given by[Sj 



r = ^v-^eM-s^d) - s c M) 



(1.1) 



to one-loop accuracy. Here S(<p c i) is the instanton action. The coefficient T> 
represents the effect of quantum fluctuations around the instanton configura- 
tion and arises from the Gaussian approximation to the functional integral. 
This is the object whose computation we will consider here. It is given in 
general form by 



the second equation relating it to the one-loop effective action. The oper- 
ators M. are the fluctuation operators obtained by taking the second func- 
tional derivative of the action at the instanton and vacuum background field 
configurations. The prime on the determinant implies omitting of the two 
translation zero modes. The first prefactor S((p c i)/27r takes into account the 
integration of the translation mode collective coordinates. Finally, the coun- 
terterm action S c t in the exponent will absorb the ultraviolet divergences of 
T>. One may also include a corresponding determinant for fermions, which 
for massless fermions is even known analytically [28, 29j. For finite masses 
is has been computed recently [9]. 

This paper is organized as follows: In the next section we present the ba- 
sic equations for the Abelian Higgs model and its instanton. The fluctuation 
operator and its partial wave reduction is presented in section [31 Two meth- 
ods for computing fluctuation integrals, one based on integrating Euclidean 
Green' s functions, as it was used in Ref. [19], and the Gel'fand-Yaglom 
method are introduced in section H] and compared. This includes the appli- 
cation for a single channel problem, as present here in the Faddeev-Popov 
sector, and to a coupled-channel problem, as present here for the gauge-Higgs 
sector. In section [5] we adress some specific problems: in subsection 15.11 we 
discuss the s-wave problem which arises do to the topological nature of the 
background field and which constitutes the main purpose of this manuscript, 
in subsection 15.21 the zero mode problem which has well-known solutions for 
both the Green's function and the Gel'fand-Yaglom method, and in subsec- 
tion 15.31 the renormalization, which here amounts to as simple subtraction. 
The numerical results, in particular a comparison of both methods and the 
final results for the effective action are presented in section [6j A summary 
and conclusions are presented in section [7J 
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2 The Abelian Higgs model and its instanton 



The Abelian Higgs model in (1+1) dimensions is denned by the Lagrange 
density (written in the Euclidean form relevant here) 

C = \{F, V ) 2 + \\D^\ 2 + \{\^\ 2 -v 2 ) 2 . (2.1) 

Here 

Ffiv = d^A u — d u A^ 

The particle spectrum consists of Higgs bosons of mass m 2 H = 2\v 2 and 
vector bosons of mass = g 2 v 2 . Usually the Higgs and gauge sector are 
coupled to a fermionic sector and displays fermion number violation violated 
by instantons. We here omit this aspect entirely. 

The model has instanton solutions which change the topological charge 

q = JL f <p X£ . (2.2) 

2n J 

If the density of instantons is sufficiently small they can be treated in the di- 
lute gas approximation and be described as separate objects with topological 
charge by q — ±1. 

A structure which exhibits such a topological charge and satisfies the 
Euclidean equations of motion is given by the Nielsen-Olesen vortex [26J. 
The spherically symmetric ansatz for this solution is given by 

Aftx) = £ -^A(r) , (2.3) 

cf) c \ x ) = v f{r)e iLp{x) . (2.4) 

In order to have a purely real Higgs field one performs a gauge transfor- 
mation 

Aft -> A^-d^ip/g 

(2.5) 
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to obtain the instanton fields in the singular gauge 



Afa) = ^(A(r) + 1) , (2.6) 

<p c \x) = vf(r). (2.7) 

In this gauge the fields take to their vacuum values as r = |x| — > oo. 
With this ansatz the Euclidean action takes the form 

oo 



rrriur V dr I \ ar I r 



+ 



rm H 



(/ 2 W-1) 2 ) ■ (2.8) 



For the case Mjj = Mw an exact solution to the variational equation is 
known [27], for which the classical action takes the value S c i = nv 2 . We will 
consider here the general case, however, for which the classical equations of 
motion 



have to be solved numerically. 

Imposing the boundary conditions on the profile functions 

A(r) const ■ r 2 4 , A(r) ^> —1 

f(r) const ■ r , f(r) '—^ 1 



(2.11) 



the Chern-Simons number is 1 and the action is finite. 

Since we will consider fluctuations around these solutions later on, a good 
numerical accuracy for the profile functions f(r) and A(r) is required. We 
have found that the method used previously by Bais and Primack [30] in 
order to obtain precise profiles for the 't Hooft-Polyakov monopole is very 
suitable also in this context. The method is outlined in Appendix lAl 
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3 Fluctuation operator and mode functions 



The fluctuation operator is defined in general form as 

5 2 S 



where ipi denotes the fluctuating fields and ipf the "classical" background 
field configuration; here these will be the instanton and the vacuum config- 
urations. If the fields are expanded around the background configuration as 
tfi = i'f + 4>i ari d if the Lagrange density is expanded accordingly, then the 
fluctuation operator is related to the second order Lagrange density via 

C n = IftMyh ■ (3-2) 

In terms of the fluctuation operators M. on the instanton and Af° on the 
vacuum backgrounds, the effective action is defined as 

Se// = IlJ^\. (3.3) 
/; 2 \detM°j K ' 

For our specific model we expand as 

1 d 



A, = Af + a^ , (3.4) 
<P = (f) cl + V . (3.5) 



In order to eliminate the gauge degrees of freedom we introduce, as in 
Ref. [31], the background gauge function 

^(A)=^ + |((0 d )*0-0V) , (3.6) 

which leads in the Feynman background gauge to the gauge-fixing Lagrange 
density 

= \{d,%f - |cv(^0 d + 4> d d^ - <p*d»4> d - 4> d d^*) (3.7) 

4V)>-^*) 2 . 



The associated Fadeev-Popov Lagrangian becomes 

C FP = l -v\-d 2 + g 2 {<P d ) 2 )v. (3.8) 

In terms of the real components f — fi + if 2 and rj = (rji + ir] 2 )/V2 the 
second order Lagrange density now becomes (omitting the superscript from 
<p d and Afj 

(C + C GF + C FP ) n = a^(-d 2 + g 2 <p 2 )a, 

+f 1 1 -(-d 2 + 9 2 Al + X(3^-v 2 ))f 1 

+ f2 l -{-d 2 + g 2 Al + g 2 <p 2 + \{< t) 2 -v 2 ))f 2 

+ f2(gA lx d fl )f 1 + f 1 (-gA fl d^)f 2 (3.9) 
+ a^(2g 2 A^(p)f 1 + a fl (2gd^(p)f2 

{~d 2 + g 2 <p 2 ) Vi + V2^ {-d 2 + g 2 4> 2 ) V2 ■ 

Specifying now the fluctuating fields (01,02,03,04,0s) as (ai, 02, fi, f2, ^12) 
the nonvanishing components of M are 

M u = -d 2 + g 2 <p 2 " * 2 1 ~ 2 ^ 2 

M 13 = 2g 2 A x <\> 
M 23 = 1g 2 A 2 $ 

M 33 = -8 2 +g 2 A 2 +\{^ 2 -v 2 ) 
-M44 = -d'+g'Al + g^ + X^-v 2 ) 
M 55 = -8 2 +g 2 <p 2 . 

The Faddeev-Popov fluctuations, labelled with the subscript 5, represent 
an single channel system, while the gauge-Higgs fluctuations, labelled with 
subscripts 1 — 4, form a 4 x 4 coupled channel system. It is understood that 
the contribution of the Faddeev-Popov operator .M55 enters with a negative 
sign and a factor 2 into the definition of the effective action. The fluctuation 
operators for the instanton and vacuum background are now obtained by sub- 
stituting the corresponding classical fields. The vacuum fluctuation operator 
for the gauge-Higgs sector becomes a diagonal matrix of Klein-Gordon oper- 
ators with masses (M W) M W) M W) M H) M w ). It is convenient to introduce a 
potential V via 

M = M° + V. (3.10) 



.M22 


= -0 2 + g 2 


Mu 


= 2gM 


M 2i 


= 2#<9 2 


M 3i 


= 


Ma 3 


= gA^O^ 
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The potential V will be specified below after partial wave decomposition. 

The fluctuation operator Ai can be decomposed into partial waves and 
its determinant decomposes accordingly, 

+00 

lndet.M= lndetM„ . (3.11) 

n=— 00 

We introduce the following partial wave decomposition for fields 

' cosy? \ e inip . ( v ( -siny? ^ pintp 



imp 



3 = E W( sine, )^T + ^( r 

n=— 00 

+00 

T MO 
^ V V2tt 

n=— 00 v 

+°° _ p in<p 

E M r ) 

=—00 

+00 

E ^( r ) 



V?2 = 



/2tt 

n=— 00 

f x D inip 



Vl2 = 



After inserting these expressions into the Lagrange density and using the 
reality conditions for the fields one finds that the following combinations are 
real relative to each other and make the fluctuation operators symmetric: 





= \{b n {r) + c n {r)) 


F?(r) 


= \{K{r) -c n (r)) 




= K{r) , 




= ih n {r) , 




= 9n(r) ■ 



Writing the partial fluctuation operators - omitting the index n in the fol- 
lowing - as 

M = M° + V, (3.12) 
the free operators M° become diagonal matrices with elements 



where (rij) = (n — 1, n + 1, n, n, n) and (Mj) = (M w , M w , M w , M H , M w ). 
The potential V takes the elements 







-1) 


v 12 




= V2m w f 




V 14 


v 22 


= v™ 




v 23 


V 24 


= -V? 4 




v 33 






n 


V 4 ™ 4 


v 55 


= ^(/ 2 - 


-1) 


v t5 







V2m w f 



A + l 



V 
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(A + l) 5 



+ 



m 



H 



(f - 1) + «4 (/ 2 - i) 



= 0. 



Choosing the dimensionless variable M^r one realizes that the fluctuation 
operator depends only on the ratio £ = Mh/Mw, up to an overall factor Myy 
which cancels in the ratio with the free operator. 

We will need in the the Euclidean fluctuation modes f"f(f, h ' 2 )i which we 
will denote as "mode functions" in the following. They satisfy 



[(M£ + ^ + Vy/^ = 



(3.14) 



The superscript a labels 4 linearly independent solutions, the subscript i 
labels the 4 components, n refers to the partial wave. The superscript + 
denotes a solution regular (i.e. exponentially decreasing) as r — > oo, the 
superscript — denotes a solution regular at r = 0. The corresponding free 
(V = 0) solutions are the modified Bessel functions B^.(Kir) with 



where Ki = \Jv 2 + Mf. 

It is convenient to rewrite the mode functions as 

The functions h^(r, u 2 ) then satisfy 



(3.15) 
(3.16) 



(3.17) 



+ 



r + B±(Kj,r) 



BtM 



[s- + K 



(3.18) 
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We note that the functions collect terms of first and higher order V n . 
The first order can be obtained by solving 



dr' 



+ 



1 BniKiT) 



^(l)a±/„ ,,2, 



(3.19) 



4 Two approaches to computing the logarithm 
of the fluctuation determinant 

4.1 Method I: Integration of the Green's function 

We consider a partial differential operator in two dimensions M. of the form 

M = -A 2 + m 2 + V(r) . (4.1) 

In the Abelian Higgs model such an operator appears in the Faddeev-Popov 
sector. The gauge-Higgs sector will be considered later. The potential V 
has been assumed to be spherially symmetric (r = |x|), so that the space of 
eigenfunctions is separable into partial waves of the form 

V> 7 (x) = it r (r)e inv (4.2) 

with 

R,(r) = -uj'R^r) , (4.3) 

where we have introduced the partial wave reduction M n of M. 0. We further 
assume that V(r) — > as r — » oo sufficiently fast, e.g. to be of finite range, 
and is nonsingular. The index 7 = (j n ,n) is a multi-index, consisting of a 
radial quantum number j n and an angular momentum ( "magnetic" ) quantum 

2 We use discrete notation for the eigenvalue spectrum, either this can be realized by 
introducing a boundary at some large value R of r or it may simply be considered as a 
formal way of presentation. 



M n i? 7 (r) 



d 1 Id n 2 , . 

\ h V(r) 

dr 2 r dr r 2 
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number n. We define the Euclidean Green's function of this operator via 

7 7 n j n ' 

= ^e^-^G n (r,r',i/ 2 ) (4.4) 

n 

satisfying 

( ^ + ^ 2 )£(x, x', z/ 2 ) = 5 2 (x - x') (4.5) 
The logarithm of the fluctuation determinant is defined as 

2 , det(/W +Iz/ 2 ) 

^ ) = ln d<AUi7) ■ (46) 

where I is the unit operator. Using the identity In det = tr In the logarithm 
can be calculated as 

j K,o) 

where u 2 fi denotes the eigenvalues of the "free" operator 

Mo = -A 2 + m 2 . (4.8) 

In order to obtain J{y 2 } we begin by integrating the Green's function 
C?(x, x, v 2 ) over x: 



(4.9) 



Next we integrate with respect to v 2 from to A 2 : 

-A 2 



^ 2 
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We now subtract the equivalent expression for the free operator Ai with the 
Green's function Q to obtain 




dv 2 [Q{v 2 ) - g (v 2 )} 




(w 7 ,o) 2 + v 2 _ 



1 



(4.12) 



n 



J(0,A 2 ) 



We now cannot longer avoid discussing the existence of the formal expressions 
we have written down. Indeed already the expression C?(x, x, u 2 ) does not 
exist. An expression that does exist is the integral over r over the partial 
wave Green's function: 



we may either sum G n (u 2 ) over n and then integrate over z/ 2 , or we may 
integrate over v 2 first and then do the sum over n. As long as A 2 is finite 
both ways lead to the same finite result. But ultimately one wants to obtain 
^7(0,00). The limit liniAz^oo J~(0, A 2 ) is naively expected to vanish, but in 
fact this is not the case, rather one finds a logarithmic dependence on the 
cutoff A 2 . It reflects the fact that the logarithm of the functional determinant 
is logarithmically divergent. This can be found by computing the leading 
Feynman graphs contributing to J{y 2 , A 2 ). 

Within each partial wave one finds, for a nonsingular potential of finite 
range, that liniAa^oo J„(0, A 2 ) is finite. But the sum J„(0, 00) is logarith- 
mically divergent, in contrast to the sum ^ n (J n (0, A 2 )). So the operations 
of summation and taking the limit do not commute. 




(4.13) 



In order to obtain 



J(0, A 2 ) = - / dv 2 [O 2 ) - G > 2 )] (4.14) 
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The logarithmic divergence can be removed by subtracting the perturba- 
tive contribution of first order in V from J„(0, A 2 ). This will be discussed in 
some detail in subsection l5.3l Then the summation over n and the integration 
over v 2 commute and the limit A 2 — ► oo can be taken. 

We finally note that for practical computations one represents the partial 
wave Green's function G n (r, r', v 2 ) by Jost functions f^(r,u 2 ) which satisfy 



(M n + v 2 )f ± (r,v 2 ) = 
the boundary conditions 



d 2 Id n 2 , . 9 
■-T-2---T- + — + V(r) + z/ 2 
ar z r dr r z 



f^u 2 ) = , 

(4.15) 



f~(r,v 2 ) oc r n for r -> (4.16) 

and 

f+(r, v 2 ) -> for r -> oo . (4.17) 
The solutions are normalized in such a way that the Wronskian is given by 

w(f + , f~) = f + (r, v 2 )^-f~(r, v 2 ) - f-(r, v 2 )-^-f + (r, v 2 ) = i . (4.18) 
These boundary conditions can be made more explicit by writing 

f~(r,u 2 ) = I n (nr)[l + h-(r^ 2 )] (4.19) 

f+(r,v 2 ) = K n {Kr) [1 + K(r,v 2 )] (4.20) 

with the condition 

lim h±(r,v 2 ) = . (4.21) 

) — >oo 

For r — > the functions h^(r, z/ 2 ) go to constants, h^(r, z/ 2 ) ~ /i^(0, z/ 2 ) + 
0(r 2 ). The Wronskian condition entails 

l + h+(0,v 2 )= 1 . (4.22) 

1 + /i n (0, z/ 2 ) 

With the mode functions defined in this way the partial wave Green's 
function is given by 

G„(r,r / )=/-(r < )/+(r > ) . (4.23) 

Having computed the solutions / for a particular set of values of n and v 2 
one know's the Green' s function for all values of r, so the r integration can 
be done to obtain G n (z/ 2 ), which is the basis for the subsequent steps. 
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4.2 Method II: the Gel'fand-Yaglom method 

The partial wave reduction of Ai 2 + v 2 is an ordinary differential operator of 
second order, and for computing the determinant of an ordinary differential 
operator there is the Gel'fand-Yaglom theorem stating that 

det(M„ + ^I) _ lim g(r,Q 



det (M ,„ + v 2 \) r->oc I 



n,0 



r, v 



from which we then obtain 



y) = l n lim ^ (r ' ^ (4.25) 



/n,' 



OV J 



r, z/^ 



Here the functions f n are identical to the functions f n of the previous sub- 
section except for the normalization. The normalization is fixed by writing 

/-(r, v 2 ) = I n {Kr) [l + h~(r, z/ 2 )] (4.26) 

with the boundary condition 

h-(0,u 2 )=0 . (4.27) 



We then have 

JJis 2 ) = ln 



1 + h n (oo, z/ : 



2\ 



(4.28) 



In the Appendix we present heuristically two proofs of the theorem, one 
along the one given in Ref. [5], and one which establishes a direct contact 
with method I. The first version of the proof is based on the condition for a 
bound state lim r ^ 00 f~(r, v 2 ) = 0. Furthermore a basic assumption is that 
Jniv 2 ) — > 1 as v 2 — > oo, i.e., that the determinant of M ra tends towards the 
one of M ni o in this limit, within each partial wave subspace. This is the case 
for potentials of finite range. But then again we have to sum over n and this 
sum will be logarithmically divergent. The renormalization for this case as 
for the Green' s function method is discussed in subsection 15.31 

We have not introduced a cutoff here and stated the theorem in its 
naive form. For the Faddeev-Popov sector this is indeed not necessary, 
the perturbative subtractions can be done in the partial waves. Alterna- 
tively, e.g. for comparing with the Green' s function method, one simply 
uses J n (0, A 2 ) = J n (0) — J n (A 2 ), where the right hand side is evaluated by 
using the naive formula. 
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4.3 Generalization to coupled channels 



Both methods can be generalized to coupled channel systems with a spherical 
symmetry. The operator M. and its partial wave reduction M now become 
N x N matrices. The solutions / ± (,rz/ 2 ) are replaced by a fundamental 
system of solutions fnt( r ^ u2 ) where index i labels the component and the 
index a labels the solution. Both indices run from 1 to N, so that the 
fundamental system can be considered to form an N x N matrix labelled by 
a and i. We have already introduced these functions and their differential 
equation at the end of section [3j 

If the Wronskian matrix within an angular momentum subspace is given 

by 



w 



a/3 



E 



a/3 



then the Green's function is given by 



,2\ tP- 



(4.29) 



(4.30) 



From this Green' s function we are again able to compute the fluctuation 
determinant using method I. We have 



GJu 2 



J rdrG n: u{r, r, z/ 2 



(4.31) 



and the subsequent steps are performed as described above. In terms of the 
functions h"^, see Eq. 13.171 we have 



G n Jr, r, v 2 ) 



KM 2 ) + K^S) 



I ni {K ir )K ni {Kir) .(4.32) 



As to the Gel'fand-Yaglom method we note that the condition for a bound 
state of the coupled-channel system is given by 



lim det r-(r,z/ 2 ) = 



(4.33) 



where the determinant refers to the matrix of the fundamental system la- 
belled by the indices a and i. Using a proof analogous to the one given in 
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Ref. [5] the theorem takes the form 



J n (z/ 2 ) = In lim i n ' l _ V - ; (4.34) 

™ det/ Q nj4 (r, i/2) 

where the matrix of free solutions can be taken as diag{J ni (/«jr)} with 
Ki = a/ v + mf. Writing the solutions f"J (r, z/ 2 ) in the form 

tij(r, v 2 ) = (fi? + Kj{r, z, 2 )) I ni (K ir ) (4.35) 

the boundary condition for the functions h"~(r, z/ 2 ) is 

lim<T(r,z/ 2 ) = (4.36) 

and the Gel'fand-Yaglom method yields 

J n (z/ 2 )=lndet{^ + ^-(oo,z/ 2 )} . (4.37) 

This has been used, e.g., in [T7|, for computing the fluctuation determinant 
for bubble nucleation in the SU(2) Higgs model at large temperature. 

We have again stated the theorem in its naive form, without cutoff. If the 
potential V n is well-behaved, we again have J n (0, A 2 ) = J n (0) — J(A 2 ). The 
situation changes if we consider the fluctuation determinant for a topological 
soliton, indeed some matrix elements of the potential V n in the gauge-Higgs 
sector are singular as r — > 0. This will be discussed in subsection 15.11 



5 Specific problems 

5.1 The case of topological solitons 

If one analyzes the operator M n for the case of the Abelian instanton one finds 
that the centrifugal barriers are modified in relation to the nontrivial winding 
number. Indeed the potential V n contains terms proportional to 1/r 2 in the 
33, 44 and 34 and 43 components. While in the topologically trivial vacuum 
sector the pattern of centrifugal barriers near r = for the gauge-Higgs 
system is {rij} = {n— 1, n+1, n, n} within the subspace of angular momentum 
n, in the instanton sector this becomes {fij} = {n— 1, n+1, n + l,n— 1}. For 
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n = we have {rii} = {1,1, 0, 0} which is distorted to {hi} = {1, 1, 1, 1}. As 
a consequence the partial waves near r = no longer behave as the Bessel 
functions J ni (/tr), and the functions h^(r, is 2 ) no longer become constant 
near r = 0. For n ^ two mixtures of the components with i — 3, 4 behave 
as 1/r and r, respectively; for n = both Higgs components behave as r. 

In Refs. [12] and EI] the method I was used in order to compute 
the fluctuation determinant for the Abelian instanton in 1 + 1 and the of 
sphaleron in 3 dimensions (the high-temperature limit of the 3 + 1 dimen- 
sional theory), respectively. The summation over partial waves was done 
before the integration over v 2 and the cutoff was sent to infinity after suit- 
able subtractions. However, if one does the v 2 integration first then one finds 
that the n = contribution becomes infinite as A 2 — > oo, even after the suit- 
able perturbative subtractions. For the Gel'fand-Yaglom method there is no 
integration over is 2 , one cannot interchange it with the summation over n, 
and one gets the partial wave contributions at once. Indeed one finds that 
the s-wave contribution, if computed naively, is ill-defined. 

The clue to the problem lies at first in the form in which the theorem 
is applied. The naive expression holds for "well-behaved" potentials. If we 
consider the version of the theorem as derived in Appendix [C] directly from 
the Green's function method, it takes the form 



det 



l + h-(oo,0) 



J n (0, A 2 ) = In =t = ^ (5.1 



det 



I + h"(oo,A 2 ) 



for a coupled channel problem. 

This form is suitable for all partial waves except for the s-wave. We had 



I + h„ (r, v 



2\ 



goes to unity. 



required that at r = the determinant det 

Relative to the free solutions, the Bessel functions, the distorted centrifugal 
barriers lead, for n ^ 0, to additional factors r and 1/r for mixtures of the 
Higgs field components 3 and 4. In the determinant these factors compensate 
and so the normalization can be maintained. However, for n = where 
centrifugal barriers change from {n,} = {1,1, 0, 0} to {hi} = {1, 1, 1, 1} two 
additional factors r appear in the Higgs sector, and the determinant behaves 
as r 2 for r — > 0. So we have a problem with the s-wave. In the Green's 
function approach the s-wave contribution is found to be finite as long as 
A 2 remains finite. The solution of the puzzle is to note that in a general 
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normalization the theorem takes the form 



det 



J n (0,A 2 ) = ln 



l + h-(oo,0) det I + h~(0,A 



det 



I + h-(oo,A 2 ) det l + h"(0,0) 



(5.2) 



Here we have set r = 
have taken the limit r 



in the additional factors, but in fact we should 
— > of their ratio. Indeed, in the limit r — > both 



det 



I + h~(r, A 2 



and det 



behave as r 2 , the ratio does have 



l + h-(r,0) 

a finite limit. So the introduction of the cutoff proves to be crucial: it allows 
for the cancellation of the factor r 2 . For suitable initial conditions at r = 
for the fundamental system h®~(r, v 2 ) the ratio can be made equal to unity 
independent of A 2 . By this modification of the numerical procedure the s- 
wave contribution can be computed, for finite A 2 , using the Gel'fand-Yaglom 
method. Once we have introduced the cutoff in the s wave we are forced 
to work with the cutoff in all partial waves, because we can let it tend to 
infinity only after summation over all partial waves. 

Even though now we have at least defined the s-wave contribution there 
remains the fact that even after subtraction of the perturbative first order 
contribution the s-wave contribution gets infinite as A 2 — > oo. This singular 
behaviour is cancelled by the sum over the higher partial waves. But this only 
works if the summation over n is performed before taking the limit A 2 — > oo. 
This will be manifest in the numerical results and will be displayed in section 
El 



5.2 Zero modes 

The m = 1 partial wave has a bound state at v 2 = 0. As the partial wave is 
degenerate with m = — 1 the bound state has a twofold degeneracy. These 
two zero modes correspond to the translational collective degrees of freedom. 
The logarithm of their eigenvalues would be infinite and has to be removed. 

For the Gel'fand-Yaglom method two methods have been proposed, by 
V. Kiselev and J.B. [16] and by Dunne and Min [TO] . The second one is 
semianalytic but has not yet been adapted to a coupled-channel problem, so 
we will use the first one: Due to the zero mode the ratio of determinants 

det(Mi + \v 2 ) 

det(M 0il + L/ 2 ) 1 ' ' 
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has a zero at u 2 = 0. The factor v 2 can be removed by taking the derivative 
of this expression at v 2 = 0. This derivative can be computed numerically 
by computing the determinant at two sufficiently small values of v 2 . So we 
have 

" det' Mi + \v 2 ~ 



. . , det Mi , d 
J JO) = In - ^ = In — 

v ; det M i dv 2 



det M i + \v 2 



(5.4) 



For the n = 1 partial wave, as for all higher partial waves, J n (0,A 2 ) is 
computed simply as the difference between J„(0) and J n (A 2 ), evaluated with 
the naive formula. For the n — 1 partial wave it is worth mentioning that 
the procedure for removing the zero mode is applied to J(0) only. In taking 
the derivative we introduce a scale, and indeed the procedure of removing 
the mode gives a dimension to the effective action, this is discussed below. 

In the Green's function method one directly computes the logarithm of 
the determinant and we have to use another procedure for removing the zero 
mode. Extending the integral — J dv 2 G\{v 2 ) from v 2 = e 2 to v 2 = A 2 the 
zero mode manifests itself by divergence lne 2 as e 2 — > 0. It is this divergence 
which has to be removed. A straightforward idea is to simply subtract the 
pole and to compute 

- £ (g^ 2 ) - 1 - G 0il ^ 2 )) + In A 2 . (5.5) 

Numerically the subtraction as we have just defined it is not suitable. The 
integral over Gi(z/ 2 ) — G ,i(^ 2 ) is convergent as A 2 — > oo, but not so the one 
over the pole term which we have subtracted. A rather simple solution of 
this problem consists in subtracting 1/u 2 only in the interval [0, 1]. Then the 
integral we have subtracted just produces the In e 2 which is to be removed 
and In A 2 is replaced by In 1 = 0. 

We note that this procedure introduces the logarithm of a dimension- 
ful quantity, while before the argument of the logarithm was dimensionless. 
Here all numerical computations are performed setting m 2 ^ = 1. So, as the 
exponential of the —1/2 log det appears in the transition rate, and as the 
translation mode is twofold degenerate, the rate will be in units M^. 



5.3 Renormalization 

As we have already mentioned the limit A 2 — > oo of J~(0, A 2 ) does not exist. 
This is due to the divergences of quantum field theory. In the present case of 
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a two-dimensional model the divergences are just tadpole diagrams propor- 
tional to 4> 2 and which have to be subtracted. We describe the procedure 
for the gauge-Higgs sector, the application to the Faddeev-Popov sector is 
obvious. 

In the Green' s function method the expression for J7(0, A 2 ), Eq. (14.141) 
is replaced by 



J(0, A 2 ) sub = - [ du 2 J2 [G n (v 2 ) - G > 2 ) - G«i 
where of course 



v 



(5.6) 



G« = J drrG n , ti (r } r } v 2 ) . (5.7) 
The diagonal components of the first order Green's function are given by 

Gn,u(r,r) = [tfj l+ (r^ 2 ) + tf,f"(r^ 2 )] I ni ( Ki r) K ni ( Ki r) (5.8) 

and the first order part of the functions h % ni is obtained by solving Eq. 
(13. 19[) which for the diagonal elements reduces to 

d 2 



dr* 



+ 



1 B^'Amr) 

-+ 2 ^ o±/ t 

r B±{Kir) 




(l)i± 



(r,i^) = VS(r) 



(5.9) 



The boundary condition is h^^(oo, v 2 ) — > 0. There is one important point, 
however: there cannot be divergences with external gauge legs, i.e. propor- 
tional to a^A^, as in a gauge theory the is no mass counter term for the 
gauge field. Indeed these tadpole contributions are compensated by second 
order terms, in scalar QED in 4 dimensions this corresponds the cancellation 
of quadratic divergences proportional to A^A^. So when computing h^^ 
using Eq. 115. 9\) the terms (A+ l) 2 /r 2 have to be omitted in Vg 3 and V^, all 
divergences are proportional to f 2 {r) — 1. 

For the Gel'fand-Yaglom method the procedure consists again in removing 
the first order part from 



Jn(0,A 2 



det 



In 



l + h-(oo,0) 



det 



I + h-(oo,A 2 ) 



(5.10) 
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The functions h^^ir,!/ 2 ) are again solutions of Eq. (15. 9p . but with the 

boundary condition h^\ l± (r, v 2 ) — > as r —>■ 0. Again the gauge field terms 
have to be omitted in the potential. The subtracted expression is simply 

l + h-(oo,0) 



det 

[U0,A 2 )] sub = In 



det [l + h~(oo,A 2 ) 



6 Numerical results 

The classical profiles were computed using the Bais-Primack method which 
is described in Appendix [S] We have used 2000 grid point for x = mwf 
in the interval [0,30], the grid was not equidistant, but the interval length 
was chosen to increase by a factor 1.005 between two neighbouring intervals, 
so as to have small intervalls at small r and larger ones in the asymptotic 
region. 

The methods of computing the fluctuation determinants has been de- 
scribed in the previous sections, this discussion already incorporates the 
numerical procedure. We have compared the two methods for computing 
functional determinants analytically. Of course this should reflect itself in 
the numerical computations. The quantity to be computed and compared is 
J„(0, A 2 ) and its partial wave sum J7"(0, A 2 ). In the following we just display 
results for n > 0, of course those for — n are identical to those for n. 

In Figured] we display the contributions of various partial waves J n (0, A 2 ) 
for small A 2 . The results of the Green's function method and of the Gel'fand- 
Yaglom method agree within drawing accuracy, and in fact to better than 
1% for all values. A difference is found for the n = 1 partial wave, due to our 
translation mode subtraction; the result are expected to agree for A 2 > 1, 
and they do. The singularity at small A 2 found in the Gel'fand-Yaglom 
results (diamonds) is due to the fact that the translation mode is removed 
from J„(0) but not from J n (A 2 ); this is correct. 

In Figure [2] we again display the quantity J n (0,A 2 ), this time for large 
A. In these results, as in the ones of the previous figure the first order per- 
turbative contribution is not yet subtracted. So neither of the contributions 
is expected to have a finite limit as A 2 — > oo, they should behave as In A 2 . 
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I 



J (0,A 2 ) 








Figure 1: Dependence of the partial wave functional determinant on the 
cutoff: small A; we display the quantity J n (0,A 2 ) as a function of A 2 . The 
solid lines are the results of the Green' s function approach, the symbols are 
those of the Gel'fand-Yaglom method; squares: n = 0; diamonds: n — 1, 
circles: n = 2; triangles: n — 5, crosses: n = 10. 



21 




Figure 2: Dependence of the partial wave functional determinant on the 
cutoff: large A; we display the quantity J„(0,A 2 ) as a function of A 2 . The 
solid lines are the results of the Green' s function approach, the symbols are 
those of the Gel'fand-Yaglom method; squares: n = 0; diamonds: n — 1, 
circles: n = 2; triangles: n = 5, crosses: n = 10. 

One sees that the s wave, n = contribution changes sign and evolves in the 
positive direction, opposite to the other ones. 

If one subtracts the first order perturbative contribution the picture changes, 
as displayed in Figure [3j Now all partial waves with n^O have finite limits 
as A 2 — > oo, but not so the one with n = 0. This is the manifestation of the 
s wave problem. At finite A 2 the sum over partial waves is convergent, and 
in the limit A 2 — > oo the singular behaviour of the s-wave is compensated by 
the other partial waves. 

Having computed the quantities J„(0, A 2 ) su b we have to do the sum over 
partial waves and then let A 2 — > oo. We find that the terms J„(0, A 2 ) su & 
behave as n~ 3 for the gauge-Higgs system and as n -5 for the Faddeev- Popov 
system. If the terms have been computed up to some n, we extrapolate by 
fitting the terms n — 5 to n to An^^ + Bn^^ + Cn' 5 and An~ 5 + Bn~ 6 + Cn~ 6 , 
respectively. Then we append the sum from n + 1 to oo using this fit. This 
procedure has been used in previous publications, it has been checked here by 
varying rh between 20 and 35 to give reliable results. The sums up to n = 35 
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50 



100 



Figure 3: Dependence of the partial wave functional determinant on the 
cutoff: large A; we display the quantity J n (0,A 2 ) as a function of A 2 after 
subtraction of the first perturbative order. The lines and symbols refer to 
the same partial waves in the previous figure. 

and the extrapolated sums are plotted in Fig. Hfor 10 < A 2 < 100. Here the 
Faddeev-Popov contributions are included. Obviously neither of these are 
independent of A 2 . The sum with the fixed upper limit n = 35 first decreases 
and then starts to increase. This is a consequence of the cancellation of 
the s-wave divergence by the other partial waves. As A 2 increases, more 
and more partial waves are necessary for this compensation, so with a fixed 
number of partial waves this cannot work. The extrapolated sum is not 
constant either, but it can be fitted to a behaviour S^ + c/ K 2 . Subasymptotic 
corrections of order 1/A 2 are expected. We consider the number as the 
asymptotic value, to be identified with the effective action. It is obvious that 
the cancellation between the s-wave and the higher partial waves becomes 
more an more delicate if A 2 increases, so it is not suitable to choose even 
higher values of A 2 to get a better estimate for S^. 
The results for the effective action 

S eff = lim ij„(0,A 2 ) su6 (6.1) 
found by the procedure we have described above are displayed in Fig. We 
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Figure 4: Cutoff dependence of partial wave sums: the sums 
l/2S n J„(0, A 2 ) su b as a function of A 2 ; circles: summation up to n = 35, 
squares: extrapolated sums; solid line: the fit —.353 + 1.037/A 2 , dashed line: 
asymptotic value —.353. The sums include the Faddeev-Popov contributions. 

also plot the results obtained in our previous publication with Torsten Daiber 
|19j . These results are consistent with the present ones within the error of 

0. 07 estimated in Ref. [19], except for the value at £ = .4. We do not intend 
to give a detailed error estimate here. From varying the maximal value of the 
angular momentum and the cutoffs we think that the error is around 0.01, 

1. e. within "drawing accuracy". The Green's function and Gel'fand-Yaglom 
methods produce consistent results within an error margin of better than 
1%. However we still have to rely on extrapolations, and this produces some 
uncertainty, which may be systematic. 

An interesting problem appears when £ = Mh/Mw > 2. The nondiagonal 
terms V34 = V43 behave as exp(— M w r). So if /^~(r, z/ 2 ) ~ exp(—K W r) it 
contributes with a behaviour exp(— (k w + M w )r to f"~(r, v 2 ). Now f£~ 
is supposed to behave as exp(— Kht). The cross term dominates over this 
behaviour for v 2 < M#(M# — 4M 2 ^)/4M^ / . This defines a real interval for 
v if Mjj > 2M\y- It found indeed that for M# > 2Mw and small v 2 some of 
the functions hf~ increase exponentially with the expected behaviour. One 
finds (numerically) that these contributions cancel in the determinant, both 
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Figure 5: The effective action; squares: our results; circles: the results of 
Ref. CHI. 



methods still produce consistent results up to £ = Mh/Mw — 2.5. However, 
this cancellation is delicate numerically and for larger values £ the numerical 
procedure breaks down, the results become inconsistent. So one would have 
to find a suitable modification of the numerical procedure in order to maintain 
numerical reliability, here we limit ourselves to the range £ < 2.5. 



7 Summary 

In this work we have addressed a problem that arises in the fluctuation op- 
erator and in functional determinant for external field configurations with 
nontrivial winding number. A modification of the centrifugal barrier factors 
nf/r 2 by the singular external field configuration necessitates modifications in 
the computation of functional determinants. While these are relatively triv- 
ial when the computation is carried out using the Green's function method, 
for the Gel'fand-Yaglom they are less obvious. We have described here both 
approaches. Indeed the handling of the problem when using the Green's func- 
tion method has led us to the solution of the problem for the Gel'fand-Yaglom 
method. It consists in introducing a cutoff, which before was unnecessary 
after suitable perturbative subtraction. 
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We have presented the numerical comparison of both methods, the results 
are found to agree within an accuracy of better than 1%. The results for the 
one-loop effective action agree with the previous calculation in Ref. [19] 
within the errors given there. 

A The Bais-Primack method 

We shortly address the problem of finding reliable classical solutions for the 
vortex background field. As we have already mentioned, we use a method 
developed by Bais and Primack [30], which has to be adapted to the vortex 
system. 

We introduce two functions F A (x) and Ff(x) with x = M w r via 



The boundary condition for these new functions are that f A (x) — > and 
Ff(x) — > as x — ► oo. For x — > they have to behave as f A (x) = l/x + 0(x) 
and Ff(x) = —l/x+ 0(x). They satisfy the differential equations 



which have been written in such a way that the differential operators on the 
left hand side 



A(r) 
f(r) 



-l + xF A (x) 
1 + F f (x) . 



(A.l) 
(A.2) 




F A (Ff - 2F f ) (A.3) 
F 2 A (l + F f ) + ^Ff(3 + F f ) . (A.4) 



F' f ' + -F'-eF f 



V A 




1 



(A.5) 



(A.6) 



are of the Bessel type. Their Green' s functions are given by 



G A (x,x') 
Gf(x, x') 



hix^Ktix^ 
Jo(ex<)K (ex > ) 



(A.7) 
(A.8) 
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and we have the integral equations 

/oo 
xdxG A (x,x') [F A (x')(Ff(x') + 2F f (x')) (A.9) 

POO 

F f (x) = F° f (x)- I xdxG f (x,x')[Fl(x')(l + F f (x') (A.IO) 
Jo 

+^Fj(x')(3 + F f (x'))-V f F°(x') . 

The functions F® und F A have been introduced in order to provide the so- 
lutions with the right boundary conditions. They have to satisfy the same 
boundary conditions as the solutions we are looking for. We have chosen 

F° A (x) = K x {x)\^- (A.ll) 
1 + CL2X 

fj(s) = -J^^-^L- . (A.12) 

with suitable parameters ai, a 2 and f\ the iteration of the integral equations 
produces solutions with an accuracy of 10~ 9 after around 150 iterations. 



B Coleman's proof of the Gel'fand-Yaglom 
theorem 

In the book by Coleman [S] the Gel'fand-Yaglom theorem is stated in the 
following way: Let f(r, v 2 ) and f°(r, v 2 ) denote the solutions of 

(M + u 2 )f(r,u 2 ) = (B.l) 

and 

(M° + u 2 )f (r,u 2 ) = 0, (B.2) 

respectively, on the interval [0, 00], with regular boundary conditions at r = 
0. Let these solutions be normalized such that 

fir v 2 \ 

1^ = 1. (B.3) 

r->0 j u (r, V A ) 
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Then the following equality holds: 

det(M + ^) _ f(r,u 2 ) 

det(M° + z/ 2 ) r-oo /o( r , u 2 ) ' y ' ' 

The argument consists of two parts: 

(i) as the bound state condition for the functions f(r,u 2 ) and f°(r,u 2 ) is 
given by /(oo, -w^) = and f°(oo, — o>° ) = 0, and as the determinants can 
be written, in the basis of eigenstates, as products with factors v 2 + u 2 and 
similarly for the determinant of M° + v 2 the right hand and left hand sides 
are meromorphic functions with identical poles and zeros. 

(ii) therefore, if furthermore both sides become unity as v 2 — > oo, they are 
identical. This condition holds for a large class of potentials, in particu- 
lar nonsingular potentials of finite range. Intuitively one expects that for 
large v 2 the potential then becomes irrelevant and that the solutions f(r, v 2 ) 
and /°(r, v 2 ) become identical, the condition my be checked by perturbative 
expansion. 

Generalized to a coupled (n x n) system the theorem can be stated in the 
following way: 

Let f (r, v 2 ) and f°(r, v 2 ) denote the (n x n) matrices formed by n linearly 
independent solutions f"(r, u 2 ) and /°°(r, u 2 ) of 

(M lJ + v 2 )f«(r,v 2 ) = (B.5) 

and 

(M»+, 2 )/; (^ 2 ) = 0, (B.6) 

respectively, with regular boundary conditions at r = 0. The lower index 
denotes the n components, the different solutions are labelled by the Greek 
upper index. Let these solutions be normalized such that 

limf(r,z/ 2 )(f°(r,z/ 2 ))- 1 = 1 . (B.7) 

r— >0 

Then the following equality holds: 

det(M + Q _ detffci/ 2 ) 

det(M° + z/ 2 ) ™detf°(r,z/ 2 ) 1 ' ' 

where the determinants on the left hand side are determinants in functional 
space, those on the right hand side are ordinary determinants of the n x n 
matrices defined above. 
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The argument goes as before, one just has to replace the bound state 
condition for the one-channel problem by the condition 



lim det f (r, v 2 ) = (B.9) 

r— >oo 

as suited for a coupled-channel problem. 



C Gel'fand-Yaglom and the Green's function 

In section H] we have have introduced two numerical methods for computing 
the functional determinant of an operator Ai, or rather of operators M n , the 
reduction of the operator M. to a subspace of definite angular momentum. 
We will here connect the two methods directly, without going back to the 
eigenfunctions of this operator which are not used in either of the methods. 
The connection between the methods has been discussed in Ref. [32J, we 
adapt their aproach to the radial and coupled-channel operators which we 
have consider here. 

We go back to the diffential equation satisfied by the mode functions 
/±(r,«) 



Id d 

r dr dr 



n 



+ + V(r) + v l 



fir, * 







(C.l) 



Taking the derivative with respect to v 2 of the differental equation for f n we 
have 



1 d d 

r dr dr 



+ -j + V(r) + v l 



d 
d^ 



f-(r,u 2 ) + f-(r,u 2 )=0 (C.2) 



Multiplying with f + (r, v 2 ) and using the differential equation for we ob- 
tain 



Id d i 2\ 
-—r—f + (r,u 
r dr dr 



— f-(r v 2 \ - f + (r u 2 ) — -—r—f~(r u 2 ) 



-G n (r, r,v 



(C.3) 
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Integrating the Green's function over r we have 



rdrG n (r, r, z/ 



00 , 1 d 
rdr- — 
r dr 



r -l-fn( r , v 2 ) 
dr 



— f-(r u 2 ) - f + (r v 2 )—r—f-(r v 2 ) 
dv 2 Jn 1 ' ' Jn 1 ' ] dv 2 dr Jn 1 ' ] 



A/„-(r, - />, v 2 )^r±f;(r, v 2 ) ) 
dv A dv A dr J 



(C.4) 



In fact we have to compute the integral over G n — Go, n and with the boundary 
conditions and normalization we have introduced the contributions of the two 
Green's functions cancel each other at the upper integration limit. Near the 
lower integration limit the functions f^(r, z/ 2 ) behave as (l+/t+(r, z/ 2 ))if n (/tr) 
and (l+h~(r, z/ 2 ))/„(/tr) with Iim P _ f0 (l+/i+ (r, z/ 2 ))(l + /?,~(r, z/ 2 )) = 1. There- 
fore the parts where the derivatives d/du 2 act on the Bessel functions cancel 
with the free contribution and we remaining term is given by 



-(1 + h+(0, ^ 2 ))^(1 + K(0, v 2 )w(K n (Kr), 7 n («r)) 



dh-(0,v 2 )/dv 2 d 
" l + h-{0,u 2 ) ~ d^ 

where we have used 



In [l + h~(0, 



w(K n (z)I n (z)) = z 



I n (z)-^-K n {z) - K n (z)-^I n (z) 



(C.5) 



(C.6) 



In subsection 14 . 21 we have introduced the functions f n (r, v 2 ) which differ from 
the /~ in the normalization. Writing 



we have 



fn ( r > V 1 ) = 

l + h-(0,u 2 )- 



1 + K(r,v 2 ) 



1 + h~ (00, v 2 



(C.7) 
(C.8) 
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So we obtain 



JnKA 2 ) = 



dv 



'2 



G n (z/ 2 ) - G >' 2 ) 



'2\ 



In 



A , d 
dv 2 —: In 

j ctz/ 2 

l + /^(0,z/ 2 



l + ^(0,z/' 2 



A 2 



= -In 



1 + h n (oo,v 



'2\ 



= In 



l + ^(oo,z/ 2 ) 



(C.9) 
(CIO) 
(C.11) 
(C.12) 
(C.13) 



1 + M°°,A 2 ) 

Generalizing to coupled channels, choosing the Wronskian derterminant 



u a P = 5 a P we have 

- / rdrG njii (r,r,v 2 ) = 
Jo 



(C.14) 



dr Jn ' 1 1 ' ' 



— f a T( r v 2 ) - f a+ (r u 2 )—r—f a -(r v 2 ) 



The subsequent reasoning about the contributions of the upper and lower 
integration limit is analogous. We now have 

-(1 + Kt(0, ^ 2 ))^(l + K~(0, v 2 )w(K^r), I ni {nr)) 

= [l + h„-(0,, 2 )]- 1 ^^ = A lndet [l + h n -(0,, 2 )] (C.15) 

and by analogous steps as for the single-channel case we arrive at 

J n (0,A 2 ) = -/ dv 2 [G n , u (v 2 ) - Go,nA^)} 
Jo 

-L 



A d 

dv 2 — lndet [l + h-{0,u 2 )} 



In 



det 


I + h-(0,0) 


det 


I + h-(0,A 2 )" 



(C.16) 



31 



As we will see it is not always possible to normalize the fundamental system 
at r = 0. A more general expression, which treats the boundaries r = and 
r — > oo in a symmetrical way is given by 



det 


l + h-(oo,0) 


det 


I + h;(0,A 2 )" 


det 


I + h-(oo,A 2 ) 


det 


i + h-(o,o) 



where now the matrix h = {h^} refers to a fundamental system with an 
arbitrary normalization. 
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